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SUI’iMARY 


BIOMOD is an operational system specifically designed to make it 
convenient for unsophisticated computer users to study models of dy- 
namic systems. The system features a high degree of interaction, 
user— oriented model-definition languages, and flexible, in-depth model 
structuring. It employs a graphics console comprising a television 

screen, a data tablet, and a keyboard. 

A BIOMOD user may represent a model by drawing block diagrams 
and handprinting or typing text; he receives immediate feedback about 
the system's interpretation of his actions. Each component of a model 
block diagram may be defined by another block diagram; this facilitates 
organizing models into meaningful substructures. A user ultimately 
defines each component block by analog-computer-like elements, alge- 
braic, differential, or chemical equations, and/or Fortran statements. 
A modeler may . thus define his model in whatever terminology is mean- 
ingful to him. Displayed curves are continually and automatically 
updated during model simulation. A user may stop the simulation and 
plot different variables, change scales, or alter parameter values, 
and then either continue simulating or revise the description of his 

model . 

This report demonstrates BIOMOD by presenting a scenario of how 
a user might describe and simulate a drug-effect model, briefly de- 
scribes the system implementation, and discusses experience with users 


Annroved 
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II. A MODEL FOR EVALUATING 


DRUG ADMINISTRATION FOLICIES 


Medications and their prescribed dosages are designed to maintain 
a critical amount of drug in the blood for a specified period of time. 
The conventional method for determining optimal dosages involves numer- 
ous laboratory experiments. If the drug effects catv be modeled, how- 
ever, a more efficient method is to experiment by running computer 
simulations. 

One technique for maintaining the prescribed amount of drug Is to 
use a .capsule containing a large number of differently coated pellets. 
The pellets dissolve .at different times, so that, as drug leaves the 
blood, it is replaced by drug released by newly dissolved pellets. 
Garrett and Lambert [8] have proposed the following model to describe 
this situation. A capsule comprises a number of pellet populations 
with different mean times of release. The rate of drug release for 
each population is assumed to be normally distributed (with the same 
standard deviation for each population) about the mean time of release 
for that population. The rate of adding drug to the body is therefore 
specified by a sum of normal distributions. The transfer of drug 
through the body is described by 


drug 




where GI refers to the gastrointestinal tract, B to the blood, and U 
to the urine and other excretory parts of the body; this means that 
drug flows from GI to B at rate k GI>B and from B to U at rate k^. 
Thus, for example, 


dDg/dt = °GI ” k B,U °B 

where D fi is the amount of drug in B, and D CI the amount of drug in GI . 

Given this description together with the requisite parameter 
values, RTOMCD can be used to simulate the model. Note that the follow 
j rr . j.-.i ar. sc --iS-- onlv cue of many possible ways or reaching 

same goal, and that BIOMOD does not force the user to take action m 



When using BIOMCD, v;e communicate via a data tablet pen and A 
keyboard. Ihe pen’s location on the tablet is always indicated by. a 
dot displayed in the corresponding location on the television screen. 
BIOMOD's interpretation of user pen actions depends on where the pen 
is placed and on what is currently displayed on the screen. We may 
handprint characters in most areas. As we write, a displayed "ink" 
track appears to flow from the pen; each time we complete a character, 
its track is replaced by a stylized character. We can change a char- 
acter by writing another over it. Some symbols are used for editing; 
for example, we may use a caret to insert text, or we may scrub with 
the pen to delete text. Some areas displayed on the screen act as 
pushbuttons; if we "push" one of these (by touching the pen down), the 
system performs the indicated action. If we push a displayed arrow, 
a continuous action takes place, such as the rescaling of a set of 
curves. Some figures can be "dragged"; if we "touch" one of these and 
move the pen, the displayed figure follows the pen's motion. We may 
type (with the keyboard) in any area where writing is possible. The 
keyboard cursor may be positioned either with the pen or with keyboard 
control keys. 

j To create our model, we enter our identification, name our model 
DRUGS, then begin constructing the model. Because it has two major 
components, we first draw two rectangles; these are replaced by styl- 
ized' function boxes. We write CAPSULE in one box and names of parts 
of the body in the other box, and then draw a flowline to connect them. 
This diagram (Fig. 2) provides no't only a picture of our model, but 
also a means of defining the two components of the model separately. 

To define the capsule component, we first push the DEFN button 
on its box. The system replaces the block diagram with a list of lan- 
guages that we may choose from to define the component. The languages 
are block diagrams, mathematical equations, chemical equations, and 
Fortran statements. We choose block diagrams so that we can define 
the capsule as a set of boxes, each representing a pellet population. 
This ability to define a box by another block diagram enables us to 
organize a model as a hierarchical collection of a number of components 
at different levels. We" draw four boxes and write PILL on the top linn 



Fig. 2 The DRUGS mode! block diagram 
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oi: o a-.i Co represent tour police populations . We name a f one cion box 
in this way whenever we anticipate using the same function repeatedly. 
Because each population is defined by a normal distribution, we 
indicate that we want to define the PILL function with mathematical 
equations. B10M0D responds by displaying a form for writing algebraic 
and differential equations. We assume that the probability of a pellet 
dissolving in an interval about time t is given by the probability 
density function 

1 „ -<tV/2o 2 


o/2tt 


Using this function to approximate the drug release rate by a deter- 
ministic variable, we write 

P » 1/ (SIGMA*SQRT (2*P I) )*EXP(- (TIME-MEAN) **2/ (2*SIGMA**2) 

The system analyzes this statement and immediately responds with the 
message 


UNBALANCED PARENTHESES 

We then add a closing parenthesis to correct the statement. We also 
realize that we should parameterize the amount of drug released by 
each pellet population, so we insert DOSAGE after the equals sign and 
scrub the 1. The display now appears as in Fig. 3. 

BIQMOD has generated separate lists of the defined and undefined 
variables; TIME does not appear because it is always the simulation 
independent variable. These lists enable us to indicate which vari- 
ables have different meanings or values each time we use the function. 
We indicate that the names (and therefore the values) of PI, SIGMA, 
and DOSAGE are the same each time we use the PILL function. This is 
because PI is a constant, and because we assume that the standard de- 
viation and dosage amount are the same for each population. On the 
other hand, we indicate that MEAN may have a different value for each 
pellet population. 
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Now chat we have defined the PILL function, we are ready to u-Je 
it to define the individual populations. Wo push a button to get back 
to our diagram of the four PILL boxes, then push the DEFN button on 
one of these. Because PILL is now defined, BI0M0D displays 

P -*■ 

MEAN f 
PI «- PI 

SIGMA SIGMA 
DOSAGE *■ DOSAGE 


for us to provide the names of the output and mean of this particular 
population. We write PI next to P -+•, to name this output PI, and write 
Ml next to MEAN to name this mean time of release Ml. We similarly 
establish the correspondence between the names of variables in the 
other three - pellet populations and the names (P and MEAN) used when 
defining the function PILL. 

We can describe the flow of the drug through the body by chemical 
equations because these are mathematically equivalent to mass transport 
equations. When we push the DEFN button on the box that describes the 
body, and select chemical equations, BI0M0D presents an appropriate 
form. According to our original model description, we would like to 
write 


GI,B 


GI 


€> Dr 


B,U 


or 


GI,B 


GI «C- 


■£=» D„ 


0 . 


B,U 


B ■<- 


-O D. 


U 


where 0. indicates that there is no backward flow, 
we iinearize each eauation, and -..'rite tne rata coc 


BIOMCD requires that 
icients and aquatic, s 
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s KGIB 0 . DGI *> DB 

S KBU 0. D3 “ DU 

Here S (for slow reaction) means that BI0K0D should derive integral 
equations from our equation. Since the pellets release drug into the 
gastrointestinal tract, we also write 


G DGI / P1+P2+P3+P4 

This Statement (with G for gain) indicates that the gain of DGI, i.e., • 
.the increased rate of change of DGI due to drug entering the body from 
outside, is equal to the sum of the rates of drug release from the four 

pellet populations. 

The model is now defined except for parameter values. When v/e 
indicate that we are ready to provide these values, BI0M0D displays 
the names of model variables and parameters in two separate lists 
(Fig. 4). Names such as DGIo indicate initial values; they are derived 
from the chemical equations by BI0M0D. We enter the values given or 
implied by Garrett and Lambert. We assume that some drug is immediately 
released ‘into the gastrointestinal tract and therefore set DGIo to 5; 
the other drug amounts are initially zero. In order to minimize storage 
requirements, BIOKQD limits the number of variables whose values are 
saved during simulation and the number of parameters whose values can 
be modified during simulation. Since this .model is small, we indicate 
that we want to save values of (and possibly plot) all the variables, 
and that we might want to modify values of all parameters except PI. 

It has taken us less than half an hour to' completely describe our 
model. We now indicate chat we would like to simulate it. BIONOD 
first produces a CS>£P/360 [2] program that describes our model and pro- 
vides for graphic display of the results. CSMP, in turn, generates a 
Fortran program, which is compiled and linked with other programs re- 
quired to run the simulation. If an error is detected at one of these 





Fig. 4 — The model 
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raescages. The time required for' the translation depends on the load 
on the (multiprogrammcd) computer; generally it is about three minutes. 

Our DRUGS model translates successfully, so the form shown in Fig. 5 
is displayed on the screen. As in the other forms, software pushbuttons 
appear across the top. The central area is for selecting numerical in- 
tegration methods, modifying parameter values, examining variable names 
and values, or plotting graphs. The areas to the left and below the 
central area are for specifying the y and x axes of the graphs. 

We expect the values of our model variables to change smoothly and 
over several units to TIME, so we choose a simple integration method — 
Simpson's method with step-size = 0.1. This is a fixed step-size method, 
so the information regarding variable step-sizes disappears. Before 


studying how to use a multi— pellec capsule, we want to ensure proper 
model behavior when there is initially some drug in the gastrointest- 
inal tract, but no capsule. To eliminate the capsule drug we push the 
PARAMETERS button to display the list of modifiable parameters in the 
central area, then overwrite the value of DOSAGE, changing it to 0. 

4 

Next we display the list of plottable variables. Because we are most 
interested in the amount of drug in the gastrointestinal tract, blood, 
and urine, we drag the names DGI, DB, and DU to the y axis. We want to 
watch the model for several simulated hours, so ws change the upper range 
of TIME (in the small box at the lower right of the central area) from 1. 
to 7. We push PLOT; now we are ready to plot DGI, D3, and DU from 0. Co 
1. against TIME from 0. to 7. hours. 

We push RESTART and the simulation begins running. We see (from 
the curves) that DB, and later DU, are being generated; the "NOW X ** " 
number changes continuously to indicate the current value of simulated 
TIME. Because CGI is plotted off scale along the upper boundary, we 
touch the pen down to stop the simulation. In order to determine the 
range of DGI, we return to the display that lists the variables along 
with their current, minimum, and maximum values. The maximum value of 
DGI is 5. (its initial value), so we write 5 over the 1 that specifies 
the upper y-axis value, then redisplay the curves. The curves are now 
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Fig . 5 — T he simulation control form with Integration AAethods 














slowly. We assume, from the curve's reasonable appearance (Fig. 6), 
chat we described at least the body component of the model correctly, 
and our choice of integration method is adequate. We see from the 
curve labeled D3 that, as expected, the drug remains in the blood for 

only a short time. 

We reintroduce the capsule drug by changing DOSAGE back to 3., 
then restart the simulation and watch the curves being continuously 
updated as it runs. Once it becomes apparent that the capsule is not 
effective, i.e., that the value of D3 drops too low, we stop the simu- 
lation. Apparently (Fig. 7), the drug is not released from the pellets 
in time to replace the drug that leaves the blood. To correct this, 
we change the mean times of release from 2., 4., 6., and 8. to 1., 2., 

3., and 4. We then rerun the simulation and get much better results. 

The amount of drug in the blood should be greater than 2.5. In order 
to determine if this is achieved, we place the pen down in the central 
area to establish an x-y meter, then drag this meter to a place where 
Y (corresponding to DB-as well as DGI and DU) is equal to 2.5 (Fig. 8). 
Our choice of means was good, but they need to be adjusted to maximize 

the total duration of capsule effectiveness. . 

While changing the means for further trials we realize that rather 
than controlling four means, we would prefer to deal only with the first 
mean and the interval between mean times of release. To reformulate 
the model in this way, we return to its description and add another box 
to the definition of the capsule component. In this box .we write 

M2 = Ml + INTVAL 

and similar equations for M3 and M4. This replaces the parameters M2, 
M3, and M4 with the single parameter INTVAL , which we set to 1. and 
mark modifiable. Once this is accomplished, we retranslate the model, 
then continue to resimulate it and change parameters until we have 
established a satisfactory drug formulation and administration policy. 




Fig. 6 — The simulation run with no capsule drug 
















pig .7 — The simulation run with the parameter values 
shown in Fig. 4 











Fig. 8 — The simulation run wifh means - 1., 


2., 3., and 4 


















